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Abstract. We discuss how to extract information about the cosmological constant from the Wheeler-DeWitt equation, 
considered as an eigenvalue of a Sturm-Liouville problem in a de Sitter and Anti-de Sitter background. The equation 
is approximated to one loop with the help of a variational approach with Gaussian trial wave functionals. A canonical 
decomposition of modes is used to separate transverse-traceless tensors (graviton) from ghosts and scalar. We show that 
no ghosts appear in the final evaluation of the cosmological constant. A zeta function regularization is used to handle 
with divergences. A renormalization procedure is introduced to remove the infinities together with a renormalization group 
equation. We apply this procedure on the induced cosmological constant A and, as an alternative, on the Newton constant G. 
A brief discussion on the extension to a f(R) theory is considered. 
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The Friedmann-Robertson- Walker model of the universe, based on the Einstein's field equations gives an explanation 
of why the Universe is in an acceleration phase. However such an expansion must be supported by almost 76% of what 
is known as Dark Energy[l]. The dark component problem results from an increasing number of independent cosmo- 
logical observations, such as measurements to intermediate and high redshift supernova la (SNIa), measurements of 
the Cosmic Microwave Background (CMB) anisotropy, and the current observations of the Large-Scale Structure 
(LSS) in the universe. The simplest candidate to explain Dark Energy is based on the equation of state P = cop (where 
P and p are the pressure of the fluid and the energy density, respectively). When CO < — 1/3, we are in the Dark energy 
regime, while we have a transition to Phantom Energy when co < — 1 . In particular, the case of co = — 1 corresponding 
to a cosmological constant seems to be a good candidate for the Dark Energy problem. Globally, this is known as 
ACDM model. Nevertheless the ACDM model fails in explaining why the observed cosmological constant is so small. 
Indeed, there exist 120 order of difference between the estimated cosmological constant and observation. Basically, the 
theoretically prediction is based on the computation of Zero Point Energy (ZPE). One possibility of computing ZPE in 
the context of the cosmological constant is given by the Wheeler-DeWitt equation (WDW)[8], which is described by 



K = %%G, Gijki is the super-metric and 3 R is the scalar curvature in three dimensions. The main reason to work with a 
WDW equation becomes more transparent if we formally re-write the WDW equation as [12] 



INTRODUCTION 
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where 
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is the volume of the hypersurface £ and 



A E = (2k) G ijkl n ij n kl - y/g 3 R/ (2k) . (4) 



Eq.(2) represents the Sturm-Liouville problem associated with the cosmological constant. In this form the ratio A c / K 
represents the expectation value of Ax without matter fields. The related boundary conditions are dictated by the 
choice of the trial wave functionals which, in our case are of the Gaussian type. Different types of wave functionals 
correspond to different boundary conditions. The choice of a Gaussian wave functional is justified by the fact that we 
would like to explain the cosmological constant (A c / k) as a ZPE effect. To fix ideas, we will work with the following 
form of the de Sitter metric (dS) 

„2_ (, Adc 2 \. 2i dr 2 



ds 2 = -[l--^r 2 j dt 2 + - — j-^ + r 2 (d9 2 + sin 2 9d<j> 2 ) (5) 
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and its counterpart the Anti-de Sitter metric (AdS) 



ds 2 = -(l + dt 2 + — % + r 2 (d9 2 + sin 2 0d<j> 2 ) , (6) 

which are different expressions of the FRW background. It is interesting to observe that Eq.(2) can be extended 
to the so-called modified gravity theories. Basically, one modifies the Einstein-Hilbert action with the following 
replacement[4] 

S=^jd 4 x^R + S ma " er -> S=^jd 4 x^f(R)+S ma " er . (7) 

It is clear that other more complicated choices could be done in place of f(R)[6]. In particular, one could consider 
/ {R.R^R^ ,R a pyS Ra ^ rS i ■ ■ •) or where G is the Gauss-Bonnet invariant or any combination of these 

quantities 1 . Even if one of the prerogatives of a f(R) theory is the explanation of the cosmological constant, we 
are interested in a more general context, where a combination of A with a general / (R) theory is considered. As a first 
step, we begin to decompose the gravitational perturbation in such a way to obtain the graviton contribution enclosed 
in Eq.(2). 



EXTRACTING THE GRAVITON CONTRIBUTION 

We can gain more information if we consider gy = gy + /zy, where gy is the background metric and fry is a quantum 
fluctuation around the background. Thus Eq.(2) can be expanded in terms of fty. Since the kinetic part of Az is 
quadratic in the momenta, we only need to expand the three-scalar curvature / d 3 Xy/g 3 R up to the quadratic order. 
However, to proceed with the computation, we also need an orthogonal decomposition on the tangent space of 3- 
metric deformations [13, 14]: 



i(a + 2V.S)* y + (Z# y +A£. 
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The operator L maps into symmetric tracefree tensors 



(^)y = V& + V^--*y(V.S), 
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hj] is the traceless-transverse component of the perturbation (TT), namely g ll hfj 



■ 0, V% 



and h is the trace of 

hij. It is immediate to recognize that the trace element a — h — 2 (V • |) is gauge invariant. If we perform the same 
decomposition also on the momentum n tj , up to second order Eq.(2) becomes 
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1 For a recent riview on f(R), see Refs.[5, 6, 7], while a recent review on the problem of /(G) and / (if, G) can be found in Ref.[9, 10]. 



Concerning the measure appearing in Eq.(2), we have to note that the decomposition (8) induces the following 
transformation on the functional measure @h t 
variable is 

1 



3>hfj@^i2>(3J\, where the Jacobian related to the gauge vector 
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+ -V ! W -R'J 



(11) 



This is nothing but the famous Faddev-Popov determinant. It becomes more transparent if t, a is further decomposed 
into a transverse part <i;J with V a ^J = and a longitudinal part §J with £jj = V a yA, then J\ can be expressed by an 
upper triangular matrix for certain backgrounds (e.g. Schwarzschild in three dimensions). It is immediate to recognize 
that for an Einstein space in any dimension, cross terms vanish and J\ can be expressed by a block diagonal matrix. 
Since detAZ? = detAdetB, the functional measure S>hij factorizes into 



®hi } = (detA£) 2 (det 



- A 2 +ViR ij Vj 
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(12) 



with ^Ay) = Ag'-' — R'-i acting on transverse vectors, which is the Faddeev-Popov determinant. In writing the 

functional measure 3>hij, we have here ignored the appearance of a multiplicative anomaly [ 1 1] . Thus the inner product 
can be written as 



/ 



S>hj 



hi 
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Nevertheless, since there is no interaction between ghost fields and the other components of the perturbation at this 
level of approximation, the Jacobian appearing in the numerator and in the denominator simplify. The reason can be 
found in terms of connected and disconnected terms. The disconnected terms appear in the Faddeev-Popov determinant 
and these ones are not linked by the Gaussian integration. This means that disconnected terms in the numerator and 
the same ones appearing in the denominator cancel out. Therefore, Eq.(lO) factorizes into three pieces. The piece 
containing is the contribution of the transverse-traceless tensors (TT): essentially is the graviton contribution 
representing true physical degrees of freedom. Regarding the vector term A£, we observe that under the action of 
infinitesimal diffeomorphism generated by a vector field £,-, the components of (8) transform as follows [13] 



(14) 



The Killing vectors satisfying the condition V,-^; + V.-^j = 0, do not change hu, and thus should be excluded from the 
gauge group. All other diffeomorphisms act on /z, 7 nontrivially. We need to fix the residual gauge freedom on the vector 
4j. The simplest choice is 4, = 0. This new gauge fixing produces the same Faddeev-Popov determinant connected to 
the Jacobian J\ and therefore will not contribute to the final value. We are left with 
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Note that in the expansion of J z d 3 x^/gR to second order, a coupling term between the TT component and scalar 
one remains. However, the Gaussian integration does not allow such a mixing which has to be introduced with an 
appropriate wave functional. Extracting the TT tensor contribution from Eq.(2) approximated to second order in 
perturbation of the spatial part of the metric into a background term gij, and a perturbation hij, we get 



A ^ = 4V j/xy/iG 1 * [(2*)*"^ (x,x) ijkl + (A L )^ (x,x) iakl 



where 



(A L /z x ) - (a l /z x ) - ARtyij + iRhfi 
is the modified Lichnerowicz operator and Aiis the Lichnerowicz operator defined by 

(AiA)y - Ahij - 2R ikjl h kl +R ik h)+R }k h k i A = -V"V a . 



(16) 
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Qijki re p resents tne inverse DeWitt metric and all indices run from one to three. Note that the term ~\B!fi^-- 
disappears in four dimensions. The propagator K^- (x,x) iak[ can be represented as 
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where h\ a 



W-L, 



( x ) are the eigenfunctions of A L . T denotes a complete set of indices and X (t) are a set of variational 
parameters to be determined by the minimization of Eq.(16). The expectation value of Aj- is easily obtained by 
inserting the form of the propagator into Eq.(16) and minimizing with respect to the variational function X (t). Thus 
the total one loop energy density for TT tensors becomes 
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The above expression makes sense only for cof (t) > 0, where 0), are the eigenvalues of A L . Concerning the scalar 
contribution of Eq.(15), in Ref.[24] has been proved that the cosmological constant contribution is vanishing for a 
Schwarzschild background. If we follow the same procedure for dS and AdS metrics, we can show that the only 
consistent value is given by Ajs = Aj S = 0. In the next section, we will explicitly evaluate Eq. (20) for a specific 
background. 



ONE LOOP ENERGY REGULARIZATION AND RENORMALIZATION FOR THE 

ORDINARY f(R)=R THEORY 



the dS and AdS metric can be cast into the following form 

ds 2 = —N 2 (r {x))dt 2 + dx 2 + r 2 (x) (d9 2 + sin 2 9d(j> 2 ) , 

where 

dx = ±- 



dr 
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with 



(A L /i x ) . . can be reduced to 



dx 2 



b(r) = —r ; 



b(rj = — r . 
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with the help of Regge and Wheeler representation [15], where we have used reduced fields of the form /, (x) =Fi(x) jr 
and where we have defined two r-dependent effective masses m\ (r) and (r) 



{r = r(x)). 



. ^W = ^(i-^) + ^'W + ^W 

In order to use the WKB approximation, from Eq. (24) we can extract two r-dependent radial wave numbers 

/(/+1) 
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kf(r,l, coi.m) = mj 
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To further proceed we use the W.K.B . method used by 't Hooft in the brick wall problem[ 1 6] and we count the number 
of modes with frequency less than co„ i= 1,2. This is given approximately by 
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(CO,) = 

Jo 



Vi{l,(di)(2l + \)dl, 



(27) 



where V; (/,©,-), i = 1,2 is the number of nodes in the mode with (/,©,), such that (r = r (x)) 

v t (l,OH) = ^ [ + dxJtf(r,l,OH). (28) 

7t J— CO * 

Here it is understood that the integration with respect to x and / max is taken over those values which satisfy kf (r, /,©;)> 
0, i = 1,2. With the help of Eqs.(27,28), Eq.(20) becomes 



1 2 / 
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This is the one loop graviton contribution to the induced cosmological constant. The explicit evaluation of Eq.(29) 
gives 

A 1 2 



s*g = pi + P2 = - g * v * - m < (30) 

where we have included an additional Alt coming from the angular integration. The use of the zeta function regular- 
ization method to compute the energy densities p\ and P2 leads to 



e \nt}(r)y/e 



i-l,2 , (31) 



where we have introduced the additional mass parameter in order to restore the correct dimension for the regularized 
quantities. Such an arbitrary mass scale emerges unavoidably in any regularization scheme. The renormalization is 
performed via the absorption of the divergent part into the re-definition of a bare classical quantity. Here we have two 
possible choices: the induced cosmological constant A or the gravitational Newton constant G. In addition, we restrict 
our investigation to the case where 

mj(r)=m 2 2 (r)=m 2 (r), (32) 
because the dS and AdS backgrounds fall in this case. 



Running Cosmological Constant 

If we adopt to absorb the divergence using the cosmological constant A, we can re-define A — > Ao + A dw , where 

Ad ,v = ^£M (33) 

e32n 2 

The remaining finite value for the cosmological constant reads 

^ = (Pi(M)+p2(M))=P e //(M,>-), (34) 

where p,- (ju) has the same form of p, (e) but without the divergence. The quantity in Eq.(34) depends on the arbitrary 
mass scale jU. It is appropriate the use of the renormalization group equation to eliminate such a dependence. To this 
aim, we impose that[17] 

1 dA (u) d TT . , 
^^T l =^PeIfM. (35) 



8nG^ dpi ^dn 

Solving it we find that the renormalized constant Ao should be treated as a running one in the sense that it varies 
provided that the scale jj, is changing 

Ap(H,r) _ Ap ( /Xq , r) m 4 Q (r) ji 
8?rG 8nG I6n 2 juq" 



Substituting Eq.(36) into Eq.(34) we find 



Ao(jUo,r) 
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If we go back and look at Eq. (2), we note that what we have actually computed is the opposite of an effective potential 
(better an effective energy). Therefore, we expect to find physically acceptable solutions in proximity of the extrema. 
We find that Eq.(37) has an extremum when 



1 
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Actually Ao (jUo, f) is a maximum, corresponding to a minimum of the effective energy. The effect of the gravitational 
fluctuations is to shift the minimum of the effective energy away from the flat solution leading to an induced 
cosmological constant. Plugging Eq.(38) into Eq.(36), we find 
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which can be set to zero when 



Ao(A,0 
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when 



H = exp 
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It is clear that the case is strongly dependent on the background choice. In this work we fix our attention on dS and 
AdS metrics written in static way 



dS and AdS background 
In the case of dS and AdS spaces, the effective masses are 



m \ (j) = m\ (r) = m\ (r) = < 
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dS Case 



(41) 



k ji+^Ads re(0,+»o) AdS Case 



The effective masses have a spurious dependence on r, which can be fixed by the extremum condition described in 
Eq.(38). This is the analogue dependence of the energy momentum tensor on the scale factor a in the Friedmann- 
Robertson-Walker model. Note that from Eqs.(41), m\ (r) can never be vanishing, except for the trivial case of 
A4rfs = 0. It is interesting to evaluate the result in proximity of the cosmological throat r c = ^/3/Aj S for 



MS 



the dS solution. In this case, Eq.(38) leads to 
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and Eq. (40) becomes 
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when 



HdS = exp 



J \is 
^ 2 ds 
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For the AdS background, we take the analogue limit of the cosmological throat, namely r 



F 2 = 
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°o, then Eq.(38) leads to 
(44) 



and Eq. (40) becomes 



Ao (P-AdS^Ads) 
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= when fl Ads = exp I - I p . (45) 



Note that at this level of approximation we are unable to distinguish contributions coming from a dS or AdS 
background. On the other hand, it is interesting to observe that if 

A ( /5<A</5 and A Ads < A Ads (46) 

then 

Prfs<Po and p Arfs < p . (47) 

This means that if we start from po at the Planck scale, we can fine tune a vanishing cosmological constant for small 
energy scales. Nevertheless to obtain this behavior, we need to assume small initial background parameters in such a 
way that condition (46) is satisfied. 



Running Newton constant 

If we adopt to absorb the divergence using the Newton constant G, we have to consider the following substitution 

1 1 m 4 (r) 



G G (p) Ae47r' 



(48) 



Nevertheless, we have to say that this procedure is not immediate for the Schwarzschild metric and related generaliza- 
tions. Indeed in this case, Eq.(48) becomes 



G^G (p)H 



1 1 /3MG (p)\ 2 1 



Ae47z; 



(49) 



which means that the divergence is not removed. Therefore, it appears that this procedure is well defined only for the 
dS and AdS cases. The remaining finite value for the cosmological constant reads now 

= (Pi (AO + Pi 00) = P T e Jf (M,0 , (50) 

where p, (p) has the same form of p, (e) but without the divergence. We eliminate the dependence on the arbitrary 
mass scale p, by imposing[17] 

A d (G7, 1 (p)) d tt , 

Solving it we find that the renormalized constant Go should be treated as a running one in the sense that it varies 
provided that the scale p is changing 

Go(p)= m ^ ■ (52) 

l + ^G (po)lnA 

Even in this case, it is interesting to consider the asymptotic part of (r) for both dS and AdS metrics. It appears 
from Eq.(52) that there is a Landau pole at the scale 



^ ex p(-X^b)-^ dSCase 



(53) 



^ ex P(-ra^)=M AdS Case 



invalidating the perturbative calculation [2, 3]. Substituting Eq.(52) into Eq.(50) we find that the expression of the 
induced cosmological constant is the same as the one in Eq.(37) with the replacement 

Aq (jup,r) _^ A(r) 
%%G ^8^G (jUo)' 

showing that in this case is the Newton constant that is running. Nevertheless, a fundamental difference in rinormaliz- 
ing the Newton constant comes from the fact that we cannot find an appropriate scale where the cosmological constant 
can be very small or eventually zero. To this purpose, we try to generalize this approach including a generic / (R) 
theory. 



EXTENSION TO A GENERIC / (R) THEORY 

It is interesting to note that Eq. (2) can be generalized by replacing the scalar curvature R with a generic function of 
R. Although a f(R) theory does not need a cosmological constant, rather it should explain it, we shall consider the 
following Lagrangian density describing a generic f(R) theory of gravity 

JSf = >/=i(/(J«)-2A) ) withf^O, (55) 

where f(R) is an arbitrary smooth function of the scalar curvature and primes denote differentiation with respect to 
the scalar curvature. A cosmological term is added also in this case for the sake of generality, because in any case, 
Eq.(55) represents the most general Lagrangian to examine. Obviously /" = corresponds to GR.[20]. The semi- 
classical procedure followed in this work relies heavily on the formalism outlined in Refs.[24, 19]. The main effect of 
this replacement is that at the scale we have a shift of the old induced cosmological constant into 



Ap(Mo^) = 1 
8tzG ^/h{R) 
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where V is the volume of the system. Note that when / (R) = R, consistently it is h (R) = 1 with 



We can always choose the form of f(R) in such a way Aq {Ho,r) = 0. This implies 



(56) 



KM_ 1 1 j-^RfiR) fjR)_ (58) 



8kG ^/h{R) 16kGV Jz v ° f'(R) 
As an example we can examine the following model 2 

f(R) =AR p exp(-aR). (59) 
With this choice, the integrated extra-potential becomes 

and the function h (R) assumes the form 

h < p) = - exp ( ~ aR) — ' {p ~ aR) Z 2 an 
K) Aexp(-aR)RP- l ( P -aR) ' K ' 



2 Several models of f(R) theories are examined in Ref.[21] 



Note that the scalar curvature is four-dimensional like the argument in f(R). One can choose, for example the 
Schwarzschild background to obtain 



\ \ ^j ^ + 2aL G P = ° 



(62) 



while for the dS (AdS) case R = ±4A and for the extra-potential one gets 

1 f R( P -aR-l) _ / f dS 

vk dX ^ 8 p-aR \ { AdS • (63) 

There exists a singularity when p = ±a4A for the dS (AdS) case, while h (R) becomes 

3Aexp(-a4A)(4A)'^'(p-a4A)-2 ^ 

, , „ , ____ . } Aexp(-a4A)(4A) p -'(p-a4A) (M) 

[ ' ) 3Aexp(a4A)(-4A) p -'(p+a4A)-2 Ms - V ' 
Aexp(a4A)(-4A)''"' (p+«4A) 



From Eq.(63) we can see that 



when 



A " ( ^ r) =0, (65) 
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Since the modified cosmological constant A' Q follows an evolution equation of the form (36), it appears that 

KM_ / \ ^InA 

\2 



8tzG \ I A Ad S ] n JL 
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(67) 



and for the value of p = 1 for both the dS and AdS background, we get a vanishing induced cosmological constant 
also at the scale jj,. However, this result is valid when we assume that the main contribution of the effective masses 
is concentrated in proximity of A c is or A^. On the other hand, when we adopt to renormalize the Newton constant, 
because of Eq.(58) and Eqs.(66), the modified cosmological constant cannot be set to zero at any scale, because there 
is no a corresponding evolution equation similar to Eq.(67). 



CONCLUSIONS 

In this contribution, the effect of a ZPE on the cosmological constant has been investigated using two specific 
geometries such as dS and AdS metrics. The computation has been done by means of a variational procedure with 
a Gaussian Wave Functional which should be a good candidate for a ZPE calculation. We have found that only 
the graviton is relevant[18]. Actually, the appearance of a ghost contribution is connected with perturbations of 
the shift vectors[13]. In this work we have excluded such perturbations. As usual, in ZPE calculation we meet the 
problem of divergences which are regularized with zeta function techniques. After regularization , we have adopted to 
remove divergences by absorbing them into classical quantities: in particular the Newton constant G and the induced 
cosmological constant A. This procedure makes these constants running with the change of the scale jx appearing in 
the regularization scheme. There are two possibilities: 

a) A is running. Then we find that it is possible to find some critical values of the renormalization scale fi where A can 

be set to zero. However, these points are strongly dependent on the background choice. The situation changes a 
little when we replace R with / (R) even if the final result depends on a case to case. In the case under examination 
of the model (59) for the value of p — 1, we find a vanishing cosmological constant at any scale. 

b) G is running. For this case, the induced cosmological constant of Eq. (50) cannot be set to zero at any scale. Actually, 

the ratio A/ (8tzGq (/x)) can be vanished. Nevertheless, the point where this happens is the Landau point, that it 
means that the procedure fails for that value. For a f(R) theory, the problem has not yet examined. 



A comment concerning our one loop computation is in order. This approach is deeply different from the one loop 
computation of Refs.[22, 23], where the analysis has been done expanding directly f(R). In our case, the expansion 
involves only the three dimensional scalar curvature. Note that with the metric (21) and the effective masses (25), 
in principle, we can examine every spherically symmetric metric. Note also the absence of boundary terms in the 
evaluation of the induced cosmological constant. 
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